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ABSTRACT 



Context. Little is known about the stellar environment and the genealogy of our solar system. Short-lived radionuclides (SLRs, mean 
lifetime r shorter than 100 Myr) that were present in the solar protoplanetary disk 4.56 Gyr ago could potentially provide insight into 
that key aspect of our history, were their origin understood. 

Aims. Previous models failed to provide a reasonable explanation of the abundance of two key SLRs, 26 Al (r 2 6 = 1.1 Myr) and 60 Fe 
(r 6 o = 3.7 Myr), at the birth of the solar system by requiring unlikely astrophysical conditions. Our aim is to propose a coherent and 
generic solution based on the most recent understanding of star-forming mechanisms. 

Methods. Iron-60 in the nascent solar system is shown to have been produced by a diversity of supernovae belonging to a first 
generation of stars in a giant molecular cloud. Aluminum-26 is delivered into a dense collected shell by a single massive star wind 
belonging to a second star generation. The Sun formed in the collected shell as part of a third stellar generation. Aluminum-26 yields 
used in our calculation are based on new rotating stellar models in which 26 Al is present in stellar winds during the star main sequence 
rather than during the Wolf-Rayet phase alone. Our scenario eventually constrains the time sequence of the formation of the two 
stellar generations that just preceded the solar system formation, along with the number of stars born in these two generations. 
Results. We propose a generic explanation for the past presence of SLRs in the nascent solar system, based on a collect-injection- and- 
collapse mechanism, occurring on a diversity of spatial/temporal scales. In that model, the presence of SLRs with a diversity of mean 
lifetimes in the solar protoplanetary disk is simply the fossilized record of sequential star formation within a hierarchical interstellar 
medium (ISM). We identify the genealogy of our solar system's three star generations earlier. In particular, we show that our Sun was 
born together with a few hundred stars in a dense collected shell situated at a distance of 5-10 pc from a parent massive star having a 
mass greater than about 30 solar masses and belonging to a cluster containing ~ 1200 stars. 
Conclusions. 

Key words, stars: general - stars: evolution - stars: rotation 



1. Introduction 

Short-lived radionuclides (SLRs) are radioactive elements with 
mean lifetimes under 100 Myr that were incorporated into me- 
teorites' primitive components such as calcium- and aluminum- 
rich inclusions (CAIs) or chondrules during the earliest evolu- 
tion phases of our solar system. Understan ding their origin ha s 
long been a major goal of cosmochemistry (iRussell et al.ll200ll) . 
as it is essential for constraining the stel lar environment of the 
Sun at its birth dMeyer & Clayton! bOOOl) . as well as for estab- 
lishing a chronology of th e solar system's first million years 
dMcKeegan & Davis(l2005l) . 

The SLRs with the longest mean lifetimes (> 5 Myr, such 
as 129 I [T129 = 23.5 Myr] or 244 Pu [r 24 4 = 115 Myr]) have 
abundances compatible with that of the expected Galactic back- 
ground owing to continuous st ar formation on kp c spatial scales 
and tens of Myr timescales (iHuss et al.l 120091) . Those SLRs 
with shorter mean lifetimes (< 5 Myr) appear to be in ex- 
cess relative to that same background abundance. Among these, 
10 Be, 36 C1, and 41 Ca can be made by solar energetic parti- 
cle irradiation of the protoplanetary disk ( Go unelle et alJl2006t 
iDuprat & Tatischefil2007l ; IJacobsen et al.ll201 ll) . For 10 Be alone, 
both an interstellar origin (iDesch et al .1120041) and a solar- wind 
implantation model have also been evoked dBricker & Ca fFee 



120101) . The origin of two key SLRs, 26 Al (r 26 = 1.1 Myr) and 
60 Fe (T6o = 3.7 Myr), remains elusive. 

After decades of measurements within CAIs (the first 
solids to have formed in our protoplanetary disk), the so- 
lar syst em's initial 26 A1/ 27 A1 ratio is well establishe d at 5.3 
x 10" 5 dMacPherson et al.lll995l: IJacobsen et al.ll2008l). though 
its h o mogeneity is s ti ll subject to debate dVilleneuve et al. I 
120091: iLiu et al.l |20H iGounelle & Russelll 12005b . The situa- 
tion is a bit more complicated for 60 Fe because its record in 
CAIs is hampered by se condary processes and by nickel nu 



cleosynthetic anomalies dBirck & Lugmairl Il988l; 



2007) 



2011 



Ouitte et al 



whose origin is far from being understood (ISteele et al 
Nickel-60 excesse s attributed to the decay of 6Q Fe 
were found in chondrules dTachibana & H uss 2003|; iTelus et all 
2011 ]), which are belie ved to have formed ~1 Myr after CAIs 
dVilleneuve et al .1120091) . Based on these data, th e presently ac- 
cepted upper limit of 6Q Fe/ 56 Fe is ~ 3 x 1 0~ 7 (lpauphas et al.l 
120081: IGounelle & MeiboirfeoiOUTelus etaDl2011h . 

To calculate the 60 Fe and 26 Al concentrations at the onset of 
solar system formation, in addition to the measured ratios pre- 
sented above, we rely on the protosolar abundances given by 
Lodders (2003). With 56 Fe/ l H = 3.2 x IP" 5 , 27 Al/ l U = 3 .5 x 
10" 6 and a hydrogen mass fraction of 0.71 (lLoddersll20()3h . we 
obtain the following concentrations in the nascent solar system 
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for 26 Al and 60 Fe: C Q [ 26 A1] = 3.3 x 10" 9 M /M = 3.3 ppb 
(parts per billion) and C Q [ 60 Fe] = 4.0 x 10" 10 M /M = 0.4 
ppb. The initial 26 Al/ 60 Fe mass ratio was thus equal to 8.2. Such 
elevated concentrations of 26 Al and 60 Fe need to be explained, 
and are the subject or the present paper. 

Though Asymptotic Giant Branch (AGB) stars have been 
proposed as a possible source of SLRs dWasserburg et al] |2006: 
iTri go-Rodriguez etaD 2009: L ugaro et al.ll2012l) . massive stars 
(M > 8 Mq) are the best candidates to account for the presence 
of 26 Al and 60 Fe in the nascent solar system. This is because 
massive stars at all stages of their evolution are present in star- 
forming regions, unlike AGB stars (iKastner & Myers|["l994l) . 

The most massive stars (M > 25 Mq) burn hydrogen for 
millions of years on the main sequence (MS) before they enter 
the short-lived Wolf-Rayet (WR) phase that precedes the super- 
nova (SN) explosion. Massive stars lose their nucleo synthetic 
products to the Interstellar Medium (ISM) via strong winds 
(during the MS and WR phase) and during the SN explosion. 
Interestingly, while 60 Fe is released only during the SN explo- 
sion, 26 Al is released during the MS, the WR, a nd the SN phases 
(iLimongi & Chieffill2006l: iPalacios et al1l2005h . 

In the class i cal S N model, first proposed by 
ICameron & Truranl (Il977l) . just after the discovery of 26 Al 
(ILeeetalJ 119761). 26 Al a nd 6Q Fe were delivered tog ether by 
a single SN dBoss et aD I201QI: lOuellette et all l2QQ9h into the 
nearby solar protoplanetary disk or prestellar core. The distance, 
r, at which an SN needs to be in order to inject an SLR at the 
solar abundance into a phase (prestellar core or protoplanetary 
disk) having a linear siz e vq reads as dCameron et al] Q995; 
iGounelle & Meibomll2008l) : 

r=?J?S^ (1) 
2 V M SLR 

where 77sn is the mixing efficiency of the SN ejecta with the re- 
ceiving phase, Fsn is the SN yield of the considered SLR, Mslr 
is the solar system mass of the SLR, A the time elapsed between 
the release of the radioactive element by the source and its incor- 
poration in the receiving phase, and r the mean lifetime of the 
SLR under scrutiny. 

Because the solar system abundance of 26 Al is far better con- 
strained than that of 60 Fe, we use the former SLR (C© [ 26 Al] = 
3.3 ppb) to calculate the maximum distance at which the SN has 
to be from either the disk or the core to deliver SLRs at the solar 
abundance. Using the SN yields of massive stars with M < 60 
M© (Huss et al. 2009) calculated by Woosley & Heger (2007) 
we obtain r < 0.4 pc for a d isk of mass 0.013 M© (lHayashi et al] 
119851: lOuellette et alJ l2Q05b and size ro = 1 00 AU with an injec- 
tion efficiency of 0.7 (lOuellette etaDl2009l) . and r < 0.6 pc for a 
core of mass 2 M© and size r = 0.05 8 pc with an injection ef- 
ficiency of 0.02 dBoss & Keiser 20loh . The receiving phases pa- 
rameters (mass and size) correspond to the observed properties 
of disks and cores adopted by the tenants of the single SN model 
(lOuellette et al.ll2005l: Boss & Keisedl2010l) . whi le the injection 
efficiencies are estimated b y the same authors (lOuellette et al] 
l2009l: lB OSS & Keisedl20l0h .' As we conservatively assumed A = 
0, the calculated distances are upper limits of the maximum dis- 
tances. In other words, if we applied a decay interval of A > 1 
Myr as required by all SN models (Huss et al. 2009), it would 
imply that the receiving phases (disk or core) lie at a few tenths 
of a parsec from the SN at most to receive 26 Al at the solar abun- 
dance. 

It is very unlikely, if not impossible, however, to find a pro- 
toplanetary disk or a dense core that close to an SN. Before they 




Fig. 1. Sketch of the model described in the text (see Sect. 2). 
Star generations #1, #2, and #3 are respectively in red, blue, and 
green. In panel (b), SNe remnants are shown with a dark con- 
tour. The reddish background symbolizes 60 Fe delivered by the 
SNe belonging to the first star generation. In panel (c), the purple 
color of the shell symbolizes the combination of 26 Al delivered 
by the single massive star wind from generation #2 and the 60 Fe 
delivered by generation #1 stars. Our Sun, born in the circum- 
stellar shell together with a few hundreds fellow low-mass stars 
(generation #3), is shown in yellow. 
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explode as SNe, massive stars carve large ionized regions in 
the ISM (called HII regions) where the gas density is to o low 
and te mperature too high for star formation to take place (iBallyl 
120081) . Observations show that even around a massive star that 
still needs to evolve for 2 Myr before it ex plodes as a SN, disks 
and cores are found several parsecs away (IHartmannll2005l) . too 
far to receive 26 Al and 60 Fe at the solar abundance. In addition, 
because SNe ejecta are vastly enriched in 60 Fe relative to 26 Al 
and their respective solar abundances (Woosley & Weaver 2007), 
all models relying on SN injection lead to a 26 Al/ 60 Fe ratio that 
is far lower than the initial solar ratio o f 8.2, unless special con- 
ditions are adopted (iDesch et al.ll201ll) . 

Following the pioneering work of lArnould et al.l (1 19971) . 
WR stars winds have recently been recons idered as a specific 
source for 26 Al alone (lArnould et al.l 120061: iGaidos et alJl2009l: 
ITatischeff et all l2010h . In the model of Gaidos et al. (2009), 
which considers injection of 26 Al at the molecular cloud scale, 
very specific conditions (such as a precise timing between the 
formation of massive stars and the Sun or stellar clusters with an 
extremel y large number of stars ) are needed. In the stimulating 
model of ITatischeff et al.l (120101) . 26 Al is delivered into a bow- 
shock-produced shell by a single runaway massive star moving 
with a velocity > 20 km/s in a dense (n ~ 100 cm -3 ) star-forming 
region. At suc h a velocity (~ 20 pc/M yr), the runaway WR star 
considered by ITatischeff et al.l (120101) would escape any dense 
star-forming region with size > 40 pc within 2 Myr, preventing 
the collection of dense gas well before the entry into the WR 
phase. In addition none of these models is generic (i.e. they fail 
to offer a common explanation for 60 Fe and 26 Al), nor do they 
constrain the solar system's geneal ogy. Finally, because of the 
rarity of WR stars (ICrowtherll2007l) . such models somehow re- 
quire a special explanation for the formation of the solar system. 

The goal of the present work is to identify a coherent model 
that accounts for the presence of 26 Al and 60 Fe in the nascent 
solar system and which is in line with the most recent astronom- 
ical observations of star-forming regions (Sect. 2). The proposed 
origin for 26 Al (Sects. 2 and 3) is entirely original and relies on 
new rotating models of massive stars. The proposed explanation 
for 60 Fe (Sect. 4) is an update of the Supernova Propagation & 
Cloud Enrichment (SPACE) model elaborated by Gounelle et 
al. (2009). Combining these two results, a generic explanation 
is offered for the presence of SLRs having a diversity of mean 
lifetimes in the early solar system (Sect. 5). In Sect. 6, we dis- 
cuss our results, focusing on the newly established solar system 
genealogy. 



2. Model sketch 

We propose to explain the SLR abundances in the nascent solar 
syste m within the framework of a common picture of star forma- 
tion dHennebelle et al.ll2009l) . namely that of sequential star for- 
mation within a spatially and temporally structured giant molec- 
ular cloud (GMC). We consider the following sequence of events 
(see Fig.[T]): 

- Time is that of the formation of a first generation (#1) of 
N\ stars in a region #1 of the GMC (see panel a of Fig. []]). 

- After a few Myr, massive stars from generation #1 explode 
as SNe and start to deliver 60 Fe into a neighboring region #2. 

- Five Myr after time (see in Sect. 4.1 how this time is esti- 
mated), an 60 Fe steady- state abundance (due to the balance 
between decay and production by SNe from generation #1) 
is established in region #2. 



Main sequence (O) 



WR SN 



\ 26 AI injection 



CAI (~ 10 4 yr) 



Collect (t* ~few 10 6 yr) 



Collapse (A c ~ 10 5 yr) 



Fig. 2. Timeline for 26 Al injection in a dense collected shell 
around a massive star, starting approximately 10 Myr after the 
birth of stellar generation #1 and lasting ~ t*. The upper part of 
the figure represents the time line of the central massive star of 
panel c in Fig. 1. In the middle part, the timeline of the 26 Al pro- 
duction and injection (arrows) is shown. Aluminum-26 is present 
in the wind some 10 5 yr after the star formation and absent when 
products of helium-burning appear at the surface (WC-type star) 
since 26 Al is destroyed in helium-burning cores. It is present 
while the star explodes as an SN. Aluminum-26 injection in the 
shell lasts during the whole collection phase and is assumed to 
stop during the collapse phase lasting a time Ac ~ 10 5 yr. The 
lower part represents the timeline of the collected shell, which 
eventually collapses to produce a protoSun and a protoplanetary 
disk in which most CAIs form. 



- At t ~ 10 Myr, a second star generation (#2, containing a total 
of N2 stars) forms in region #2, partly due to the compressive 
action of the ge neration #1 SNe Shockwaves onto the molec- 
ular cloud gas (iPreibisch & Zinneckerll2007l) (see panel b of 
Fig. ID- 

- From t ~ 10 Myr and for some Myr, the wind of one or 
two massive stars from generation #2 will collect ISM gas 
to build a dense shell surround ing an HII region of radius 
5-10 pc dDeharveng et al.ll201Qh . That collected shell, which 
contains 60 Fe originating in the SNe of the first generation, 
will be wind-enriched in 26 Al via efficient turbulent mixing 
(iKoyama & Inutsu ka 2002) during a time t+ lasting a few 
Myr (see Fig. [2]). 

- At t ~ 10 Myr + t+, aluminum-26 delivery ends when 
the dense shell fragments and collapses via a diversity of 
gravito-turbulent mechanisms, such as grav itational instabil- 
ities a nd ionization of a turbulent medium dDeharveng et al.l 
120101) . to form a third generation star cluster (#3). The col- 
lapse phase lasts Ac ~ 10 5 yr. Our Sun belongs to that third 
generation of stars (see panel c of Fig. [T]). CAIs formed in the 
protoplanetary disk surrounding the protoSun on a timescale 
of a few 10 4 yr (Jacobsen et al. 2008) contains 26 Al from 
the wind of the generation #2 massive star and 60 Fe from the 
generation #1 SNe (Fig.0. 

In the following, we show that such a scenario presents many 
realizations that permit a good fit of the quantities of 26 Al and 
60 Fe at the solar system's birth. In that respect, it provides a 
more natural explanation than previous models, because it does 
not require special and/or unlikely conditions, and leads to un- 
derstanding how the solar system forms within a broader per- 
spective. The key point is that 60 Fe comes from the SNe of the 
first generation of stars, while 26 Al comes from the wind of a 
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single massive star belonging to a second star generation. This 
decoupling between the sources of 60 Fe and of 26 Al is essential 
to explain the observed 26 Al/ 60 Fe ratio when the solar system 
formed, and it does appear as a signature of the sequential star 
formation process. 



3. Origin of 26 Al in the solar system 

3.1. 26 Al yields in massive stars winds 

The 26 Al yields during the wind phase are calculated using 
new models for massive stars (Fig. [3]). Th ese models present 
two major improvements ov er previous ones (iLimongi & Chieffil 
[2006k iPalacios et al] I2005I) considered so far by modeler s 
(lArnould et al.lll997l : iGaidos et al.ll2009k iTatischeif et afll2010h . 
First, t hey are built on new ly determined initial solar abun- 
da nces (lAsplund et al.ll2009l) . The solar abundances determined 
by lAsplund et al.l (I2009I) are similar to those shown by present- 
day massive stars in the solar neighborhood, and this composi- 
tion is also quite similar to those expected in massive stars at 
the fo rmation of the solar system (see Table 5 in lAsplund et al.l 
l2QQ9h . 

Second, these models accou nt for the effects of a xial rotation 
and improved mass loss rates dMevnet et a l. 2008). Rotational 
mixing is accounted for as in Palaci oset al.l ([2005). More pre- 
cisely, rotational mixing accounts for the effects of meridional 
currents and she ar turbulence. The ve rtical shear diffusion was 
taken as given bv lTalon&Zahnl(ll997l). and the horizontal shear 
diffusion was taken as in Maede r nl2003h . A moderate overshoot- 
ing is included. The radii of the convective cores given by the 
Schwarzschild's criterion are increased by 0.1 H p , where H p is 
the pressure scale height estimated at the Schwarzschild bound- 
ary. The mass loss rates are those of Vink et al. (2000, 2001) and 
of de Jager et al. (1988) outside the application domain of the 
prescriptions by Vink et al. (2000, 2001). During the Wolf-Rayet 
phase, we use the mass loss rates of Nugis & Lamers (2000). The 
impact of rotation on the mass loss rates are accounted for as in 
Maeder & Meynet (2000). 

These models reproduce reasonably well the observed num- 
ber ratio in the solar vicinity of WR to O-type stars, of nitrogen- 
rich (WN stars) to carbon-r ich (WC) WR stars, and of the tran- 
sition WN/C to WR stars dMevnet et al.ll2008h . Because these 
ratios are sensitive to both rotation and mass loss, these good 
agreements provide some support to the way these processes are 
included in the present stellar models. 

We stress here that rotational mixing allows some products 
synthesized in the core to appear at an early stage of the evolu- 
tion of massive stars. For instance, while in non-rotating models, 
one has to wait until the H-rich envelope is nearly completely re- 
moved by the stellar winds in order for the star to eject 26 Al, in 
rotating models the surface is already 26 Al-enriched during ear- 
lier phases thanks to diffusion of 26 Al from the core to the sur- 
face through the (shear unstable) radiative layers. This lengthens 
the period during which the wind is 26 Al-enriched (Fig. [3]). This 
period is no longer reduced to the WR phase as in non-rotating 
models, but also covers part of the MS phase. 



10 ' 



3.2. 26 Al abundance in the collected shell 

As already explained above, the wind of one or two massive 
stars from generation #2 (containing a total of N2 stars) will 
collect ISM gas to build a dense shell. The dependence with 
time of aluminum-26 concentration in the shell, C[ 26 A1], reads 
as C[ 26 A1] = ^wind^wind[ 26 Al](0/M she ii, where n wind is the mix- 




Time (Myr) 



Fig. 3. Integrated 26 Al abundance produced by the stellar wind of 
massive stars as a function of time. Labels indicate stars' initial 
masses (in M©). These numbers take 26 Al radioactive decay into 
account (see Eq. 1 of Arnould et al. 2006). 
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Fig. 4. Thick lines show the acceptable model solutions, i.e. the 
values of t+ and 7/ w i n( j for which the third generation stars have a 
solar 26 Al abundance. By definition /7 w inci cannot be larger than 
1 . The filled circle shows for information the onset of the WR 
phase. The dotted line shows the 60 M© model with Ac = 0.5 
Myr instead of 0.3 Myr (see text). 



ing efficiency of the wind with the shell, M s h e ii is the mass of the 
collected shell and J P w i n( j[ 26 Al] is the quantity of non-decayed 
26 Al (i.e. survivor) present in the total wind ejecta of the mas- 
sive star at time t. Imposing that the Sun forms in that shell after 
the collect and collapse phases, i.e. at a time + Ac (Fig. [2]), we 
have 



C Q [ 26 A1] 



^wind^wind [ 26 Al] (t+) 
6 



-A C /r 2 6 



(2) 



where C©[ 26 A1] is the 26 Al abundance at the onset of solar sys- 
tem formation. The last term of the equation accounts for the 
decay of 26 Al during a time Ac, which is the duration of the final 
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Cluster size 



straints exposed above. We use the IMF of iKroupa et al.l (Il993l) 
mimicked by the generating function: 

M = 0.08 + (0.19£ L55 + 0.05£ a6 )/(l - £ a58 ), (3) 

where <f is a rando m number to be chosen between and 1 
(iKroupa et al.l Il993h . We consider only those distributions for 
which the most massive star is less massive than 150 Mq, a 
likely upper limit for stellar masses (IWeidner & Kroupall2006l) . 
Using that generating function, the fraction fsN of stars more 
massive than 8 M© (i.e. which will go SN) is 2.4 x 10" 3 , and the 
average stellar mass is M* = 0.5 M©. 

The calculated distribution probability is shown in Fig. [5] 
It is bell-shaped, and clusters with size -1200 stars have the 
highest probability (-15%) to satisfy the dual constraint evoked 
above. In 10 % of the positive runs, clusters with 1200 stars con- 
tain two stars more massive than 32 M©, which would both con- 
tribute to the 26 Al enrichment. In such occurrences, the increased 
amount of 26 Al would relax the (already loose) constraint on 

?7wind • 



Fig. 5. Probability V (in %), as a function of the cluster size N, 
to realize the double condition: 1) the number of massive stars 
(M > 8 M©) is below nsB = 5 and 2) one star at least has a mass 
over 32 M©. For each point, the IMF (see text) was simulated 
10000 times. 

collapse phase, when the 26 Al-carrying wind is conservatively 
assumed not to penetrate the shell beca use it has become to o 
dense (Fig. Et . With M, he]] = 1000 M^ (iZavagno et al.ll2007l) . 
A c = 0.3 Myr (I Ward-Thompson et al.ll2007l) and C©r 26 All = 3.3 
ppb (see Sect. 1), Eq. (2) can be solved for 7/ w i n d and t+ (Fig. [4]). 

As long as M > 32 M©, there is a wide diversity of solutions 
to our problem (Fig. |4]). An acceptable solution (among many 
others) is for example enrichment by a 85 M© star during t* = 3 
Myr with an injection efficiency 7/ w i n( j = 0.07. The acceptable 
duration of enrichment (t*, see Fig. [2]) lasts from 0.65 to 6.2 Myr, 
while 77 w i n d ranges from 2.1 x 10" 2 to 1. For many solutions, 
26 Al-enrichment ends before the onset of the WR phase; i.e., 
there are many values of t* that are lower than twR, where twR is 
the time of entry in the WR phase (see Fig.©. 

3.3. Constraining generation # 2; the parent cluster size 

Unlike other models, we can constrain the size of the parent 
cluster of our solar system. If cluster #2 con tains many massive 
stars, it is likely that a super wind bubble (iParizot et al. | l2004l) 
will form and that it will open into the ISM (iGiidel et al.ll2008l) 
provoking leakage of 26 Al atoms. Nice spheri cal shells are ac- 
tually observed around single massive stars (iDeharveng et al.l 
120101) rather than around large clusters (IGiidel et al.ll2QQ8h . We 
can therefore constrain the size (N2) of the cluster #2 in im- 
posing two requirements: 1) Cluster #2 contains fewer than 
five stars more massive than 8 M© because superbubbles form 
when the number of massiv e stars is larger than nsB = 5 
(iHigdon & Lingenfelterll2005l) . 2) Cluster #2 contains at least 
one star more massive than 32 M©, the minimum mass for which 
Eq. (2) has a solution (Fig. |4j). We call ^(N) the size-dependent 
probability distribution of a cluster containing N stars and satis- 
fying conditions 1) and 2). 

To calculate P(N) we generate, for a given cluster size, 
10000 realizations of the stellar initial mass function (IMF), 
and count the fraction of realizations which satisfy the two con- 



4. Origin of 60 Fe in the solar system 

4.1. Updating the SPACE model 

In our generalized model, 60 Fe originates in SNe from the 
first generation of stars (see panel a in Fig.l), as proposed in 
the Su pernova Propagation And Cloud Enrichment (SPACE) 
model (iGounelle et al.l I2009I) . The calculations presented here 
are modified when tak ing the new (longer) mean lifetime of 
60 Fe (iRugel et al.ll2009l) into accou nt, and when adopting the SN 
yields of IWoosley & Hegeri (I2007I) rather than a combination of 
yields as in Gounelle et al. (2009). 

We consider a set of N\ stars coevally formed at time and 
distributed in mass according to the initial mass function. The 
time-dependent amount of 60 Fe ejected by SNe in the interstellar 
medium, PsN[ 60 Fe](A/i, t), writes as 

_ {t _ H) 

Ps N [ 60 Fe](M, = J] ^s Nl ( 60 Fe) e~ *o , (4) 
i=\ 

where FsN^Fe) are the 60 Fe yields of the z th SN, t { is the explo- 
sion time of the z th SN, T60 is the 60 Fe mean lifetime, and A^sn 
is the total number of SNe (i.e. stars more massive than 8 M©) 
from ge neration #l.The mass ive stars lifetimes are those calcu- 
lated bv lSchaller et all (Il992l) . 

To fix ideas, we adopt N\ - 5000 stars and perform the cal- 
culation 1000 times to take the stochastic nature of star for- 
mation into account. The evolution of PsN[ 60 Fe](5000) with 
time is shown in Fig. [6] It takes roughly 5 Myr to establish a 
steady state abundance of 60 Fe (due to the balance of SN pro- 
duction and decay). We define the background (steady-state) 
value as P S N[ 60 Fe](M) = P S N[ 60 Fe](M, 10 Myr), i.e. at the 
onset of star formation of generation #2. In our fiducial case, 
A>N[ 60 Fe](5000) = 2.1 x 10" 4 M© (Fig. 0). For each cluster 
size, 10000 realizations of the IMF were made to check the 
dependency of PsN[ 60 Fe] with N\. It is obviously linear and 
we find P S N[ 60 Fe](M) = p SN [ 60 Fe] x N u with /?s N [ 6() Fe] = 
4.2 x 10" 8 M©. 

A GMC region having a mass M 2 receiving PsN[ 60 Fe] solar 
masses of 60 Fe will have an 60 Fe concentration of C[ 60 Fe] = ysN 
£ SN [ 60 Fe] / M 2 = rsN x /?sN[ 60 Fe] Ni/M 2 , where y SN is the dilu- 
tion factor of the ejecta in region #2 and N 2 the number of stars 
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Fig. 6. Psn[ Fe](Afi, t) for N\ = 5000 stars, averaged over 1000 
realizations of the IMF. The slight difference with Fig. 2 of 
Gounelle et al. (2009) is due to the longer mean lifetime of 60 Fe. 



formed in region #2. With M 2 = M^A^/e where e is the star for- 
mation efficiency and M* is the average mass of a star (see Sect. 
3.3), one obtains the simple relationship 



C[ 60 Fe] = rsN /?s N [ 60 Fe] N x 



(5) 



With RsN rFe 



_ = 4.2 x 10" 8 , 6 = 30 % (Elmegree3|2QQl 
lLada & Ladall2003h and M* = 0.5 M Q (see Sect. [33]), it gives 
the relationship 



C[ 60 Fe] = 25.2 10" 9 y SN ^ M /M 



(6). 



Imposing that, at the onset of star formation in generation #3, 
the 60 Fe concentration established by the N\ stars from the first 
generation of stars is equal to that of the nascent solar system 
(C©[ 60 Fe]), equation (6) reads as 



C o [ 60 Fe] = 25.2 rS N^^ /T60 ppb, 



(7) 



where t is the time elapsed between the onset of star generation 
#2 and of generation #3 (that of our Sun). 

4.2. Constraining generation # 1: the grandparent cluster 
size 

The quantity t in Eq. (7), which is the duration between the onset 
of star generation #2 and of generation #3 can be written as the 
sum of t+ and of A c . With N 2 = 1200 stars (see Sect. 3.3) and 
0.65 < U < 6.2 Myr (see Sect. 3.2), we obtain 26 < y S N N\ < 
108 or, to account for a possible revision of the initial solar sys- 
tem content of 60 Fe, expressed as a function of the initial ratio 
( 60 Fe/ 56 Fe) : 



( 6Q Fe/ 56 Fe) ( 6Q Fe/ 56 Fe) 
26 — — ttt^i < 7sn N\ < 108 „_ 1/W7 • (8) 



3x10" 



3x10" 



The range of possible solutions (N\, ysn) is shown in Fig. 
13 A first generation cluster with N\ = 5000 stars and ysN = 
0.02 is a possible reasonable solution. If, as recently suggested 



Fig. 7. The dilution factor, ysN vs. iVi, the number of stars 
in star generation #1. The allowed range of parameters is de- 
picted in white and defined by ysN < 1 (by construction) and 26 
< 7sn N\ < 108 (see Sect. 4.2). 



(iMoynier et al.ll201 ll) . the solar system initial abundance of 60 Fe 
was lower than the one usually accepted, ysN or N\ could ac- 
cordingly be decreased. 

5. The origin of SLRs: a generalized model 

It is therefore possible to identify the sequence of events that ac- 
count for the presence of SLRs with a variety of mean lifetimes 
in the nascent solar system (Fig. 1). 

- SLRs with the longest mean lifetimes (> 5 Myr) and with 
abundance compatible with a Galactic background origin, 
come from a large diversity of stars on the (kp c) scale of 
star c omplexes on a few 10s of Myr timescales dHuss et al.l 
120091) . 

- 60 Fe in the solar system was synthesized by a handful of SNe 
on the GMC (10s of pc) scale on a 5-10 Myr timescale. 

- 26 Al was finally delivered by the wind of a single massive 
(M > 32 Mq) star on a few Myr timescale and on a 5-10 pc 
spatial scale. 

The variability of temporal scales recorded by SLRs' mean 
lifetimes corresponds to the diversity of spatia l scales at which 
star formation is observed (lElmegreenl 120071) : the longer the 
mean lifetime, the larger the scale of injection in the ISM. 
This spatio/temporal variability is simply th e reflection of se - 
quential star formation in a hierarchical ISM (lElmegreenll2007l) . 
Sequential star formation is an important me chanism of star for- 
mation that has been known for a long time (lElmegreen & Ladal 
1 19771: lElmegreenl [20071: iHennebelle etaDl2009l) . It is observed 
in a fo ssilized form in star-formi ng regions such as Orion ([Ball v 
2008) or Scorpius-Centaurus dPreibisch & Zinneckerl 120071) . 
where OB subgroups with different ages are interpreted to have 
arisen from propagating star formation. It has also been caught 
in the act at the b orders of HII regions dDeharveng et al.ll2010l: 
ISnider et alJkOQg^ as well as on larger scales in the Carina star- 
forming region Ichen et alJl2007h . 

Though occurring on a wide range of spatial scales, a sim- 
ilar mechanism accounts for the origins of 26 Al and 60 Fe. For 
both radionuclides, ISM gas is collected by the same agent as 
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the one that carries the SLR (SNe ejecta or massive star wind). 
After it has become dense enough, the collected gas is unstable 
to gravity and collapses to eventually form new stars. The SLRs 
that were contained in the collecting agent end up in the new 
generation of stars. In other words, given that the SNe shocks 
and massive star winds that contribute to the building of dense 
phases of the ISM also carry radioactive elements, it is not sur- 
prising that molecular clouds - and stars formed within them - 
are radioactive. 

6. Discussion 

Our model therefore provides a natural explanation for the elu- 
sive presence of 26 Al and 60 Fe in the nascent solar system. It is 
generic bec ause, both for 26 Al and 6Q Fe, a collect-and-collapse 
mechanism (lElmegreen & Ladal 19771) is at work, though operat- 
ing on different scales. It does not suffer from the other difficul- 
ties encountered by the previous models. Because the 60 Fe and 
26 Al have a decoupled origin (SN and wind, respectively), our 
model is not faced with the pr oblem of the low 26 Al/ 6Q Fe ratio 
met b y the single SN models (lOuellette et al.ll2Q05b (Boss et al.l 
l201Qh . In addition, when the massive star goes SN, the disk is 
far enough awa y (5-10 pc) to av oid the disruption caused by the 
SN Shockwave (Chevalier 2000) or over-injection of 60 Fe, unlike 
the model of Gaidos et al.l (120091) . More importantly, it agrees 
with the astronomical observations of star-forming regions and 
the accepted astrophysical mechanisms of star formation. 

Our model leaves room for some imperfect mixing to ac- 
count for the injection of 26 Al by a stellar wind (?7 wind ^_Lb un ~ 
like previous models that as sume perfect mixing (iGaidos et al.l 
120091: iTatischeff et"aD l2010h . Though it would be desirable to 
perform numerical simulations to calculate dust and gas injec- 
tion efficiencies from the first principles, the values required in 
our model (as low as 0.02) can be reached easily. In fact, it was 
shown that the injection efficiency, at least for SNe ejecta, was 
multiplie d by a factor of 3 when density was decreased by a fac- 
tor of 2 (iBoss et al.ll2010h . By construction, the collected shell 
has a lower density than that of a coll apsing core, for which in- 
jection efficiencies are as high as 0.02 (iBoss & Keis er 2010). 

The model is quite robust to parameter changes. First of 
all, it does not pretend to give a unique solution but rather a 
range of possible solutions. This diversity of solutions is in- 
tellectually satisfying since it would be fairly vain to pretend 
to exactly identify the star(s) responsible for the solar sys- 
tem 26 Al and 60 Fe enrichment relative to the Galactic back- 
ground. It also means that there is ample room for parame- 
ter changes. Second, we have been quite conservative in the 
choice of our model inputs. This is exemplified by the adop- 
tion of new stellar models with metallicity lower than the one 
adopted by previous works. This lower metallicity implies lower 
initial content of 25 Mg, which is the nucleus whose transfor- 
mation produces 26 Al during the core H-burning phase, and 
thus our models produce smaller q uanti ties of 26 Al than the 
models of iLimongi & Chieffil (120061) and iPalacios et al.l (12005 ) 
used in previous wor k (lArnould et al .111 9971: IGaidos et alJ 2009 : 
ITatischeff et al.1120101) . On the same line, we have conservatively 
considered that, during the final phase of collapse lasting Ac 
= 0.3 Myr, the shell wa s closed to 26 Al i njection, unlike, for 
example, the models of IBoss eta! ] (120101) . Adopting an even 
longer isolation time would not change our conclusions (see the 
model with Ac = 0.5 Myr in Fig. [4]). The range of values found 
for t* is compatible with the observed timescales for the onse t 
of star formation in a collected shell (IDeharveng et al.l 20101) 
The shell mass agrees with observations (Z avagno et al. 1 12007b 
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Fig. 8. Probability P (in %), as a function of the cluster size N, 
to realize the double condition: 1) the number of stars with mass 
higher than Msn is smaller than nsB and 2) one star at least has 
a mass greater than M m i n , where the values nsB, M m i n , and Msn 
are set to vary. For each point, the IMF (see text) was simulated 
10000 times. Filled circles: same as Fig. i.e. nsB = 5, M^n = 
32 Mq and Msn = 8 M©. Gray circles: all being the same but 
Msn = 10 Mq. Green circles: all being the same but nsB = 4. 
Green squares: all the same but nsB = 6. Yellow circles: all the 
same but M m i n = 40 Mq. Yellow squares: all the same but M m i n 
= 25 M Q . 



IDeharveng et al.ll2008l) though it could be a fa ctor of two higher 
dZavagno et al.ll2007l: IDeharveng et al.1l2008l) . In any case, the 
assumed mass of the shell is 10 5 times higher than the disk of 
Ouellette et al. (2005), 500 times more than the core of Boss & 
Keiser (2010), and four times greater than the most massive bow- 
shock shell of Tatischeff et al. (2010). It means that the probabil- 
ity of enrichment (not even taking the astrophysical context into 
account, Gounelle & Meibom 2008), for a given massive star, 
is far higher in our collect-injection-and-collapse model than in 
any other model, considering that a single massive star is respon- 
sible for the solar system's 26 Al. 

On a slightly different note, our reasonin g does not critically 
depend on adopting the IWoosley & Hegerl (120071) SN yields. 
These wer e chosen becau se they are based on the solar abun- 
dances of lLoddersI (120031) . which a re almost iden t ical to the 
newly determined sol ar abundances oflAsplund et al] (120091) . un- 
like the older ones of lAnders~& Grevess ej (1 19891) used b y other 
modelers of SN yields, such as Limongi & Chieffi (2006). These 
new abundances, calculated using 3D hydrodynamical models of 
the solar atmosphere, are supported by a high degree of internal 
consistency and by agreement with values obtained in the so- 
lar neighborh ood, though some dis crepancy with heliosismology 
data persists (lAsplund et al.ll2009|). It is also i mport ant to note 
that, though uncertain, the Wooslev & Hege3 (120071) SN yields 
are probably correct within a factor of a few since the IMF inte- 
grated 26 Al/ 60 Fe ratio of SNe is comparable (within a factor of 
a fe w) to the one obser ved in the ISM by the INTEGRAL satel- 
lite ( W ang et ai1l2007l) . Finally, adopting different yields (such 
as those of Limongi & Chieffi 2006) would change the results 
by a factor of a few, not critical in the reasonings of Sects. 1 or 
4, especially given the debates concerning the initial abundance 
of 60 Fe (see below). 
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An important aspect of our model is that it constrains the 
solar system genealogy; i.e., it provides estimates of the num- 
ber of stars contained in the parent and grandparent clusters (N2 
and N\, respectively). These cluster sizes are not constrained in 
the same fashion. In either of thesecases, it is important to keep 
in mind that we onlyhave access to an order of magnitude, as 
expected for a stellar nursery that was active 4.5 Gyr ago and 
whose only memory is preserved in meteorites' most primitive 
phases. 

The best estimate of N2 depends on a variety of numerical 
values as described in Sect. 3.3, namely the minimum number 
of stars, nsB, needed to initiate a super wind bubble expansion 
(identified to the minimum number of stars needed to create 
an SN-blown superbubble since the energy contained in stellar 
winds is comparable to the energy contained in SNe explosions, 
Parizot et al. 2004) and the minimum mass for which Eq. (2) has 
a solution. It also depends on the minimum mass for which a star 
goes SN which is = 8 M© but is somehow model depen- 
dent (ISmarttll2009l) . To check that dependency, we explored the 
variations in !P(N) with these parameters (Fig. [8]). Even allow- 
ing a relatively large and unlikely exploration of the parameter 
space, it is easy to see that the most probable value of N2 always 
lies between 1000 and 2000. This exploration of parameters also 
encompasses IMF formulations other than the one in Kroupa et 
al. (1993). 

Given the dependency of the size of generation #1 on the ini- 
tial abundance of 60 Fe and on ysN (see Eq. (8)), it is obvious that 
the value of N\ is at present poorly known. It is not surprising 
that the size of the grandparent generation is less precisely con- 
strained than the parent generation. As in any family, the further 
back in time one goes, the more difficult it is to retrieve precise 
information on ancestors. More knowledge will be gained on N\ 
when a nalysts have converged on the initial sol ar system value 
of 60 Fe (lOuitte et al.ll2010l:lMovnier et al.ll201lh . along with the 
development of numerical simulations that will help to better 
constrain the value of ysN- 

The model can naturally account for some heterogeneity in 
the initial distribution of 26 Al. Injection takes place on only 
one side of the dense shell. Because some time will be needed 
for 26 Al carriers (gas or solid) to diffuse within the shell en- 
tire volume, it is expected that the firs t solids to form close 
to the protostar, as expected for CAIs (iKrot et al.1 120091) . will 
not contain 26 Al, as is observed in some meteoritic samples 
dLiu et al.ll2009tlMakide etZll201ll:IWeber et al.lll995h . Though 
the exact diffusion timescale is difficult to quantify without en- 
tering desired complex numerical simulations, it is expected to 
be close to the crossing time of the core, i.e. within the range of 
0.5-1 Myr (iBergin & T afalla 2007), mar ginally compatible with 
the co llapse timescale of a few 0.1 Myr (W ard-Thompson et al.l 
120071) . Homogenization mechanisms will take place d uring 
the di sk phase following the mechanisms presented by iBossI 
(l201ll) . Whether that homogenization goes t o comple t ion is 
at present unknown (IVilleneuve et al.l 120091: iLiu et al] l2012t 
iGounelle & Russelll l2QQ5h . Given that 60 Fe production halts 
some Myr and one star generation before the Sun's forma- 
tion, its distribution is expected to be homogeneous. At present, 
there is no experimental consensus on the homogeneous vs. 
heterogeneous distribution of 6Q Fe in the early solar system 
dDauphas et al.ll20Q8b [Quitte et al.ll201Qh . 

Our model produces no undesired collater al enrichments 
(iNichols et al.lll999b IGounelle & Meibomll2007l) . Besides 26 Al, 
massive stars winds do not c ontain other SLRs un til the very last 
3 xlO 5 yr of the star's life (lArnould et al.ll2006l) . during which 
26 Al injection is not considered. Given the high dilution factors 



of the wind needed to account for 26 Al (< 10" 3 ), the oxygen iso- 
topic composition of the shell is not modified. The shift in the 
oxygen isotopic composition is calculated to be less than a few 
per mil for all the considered models, which is far below our 
detection capability. 

7. Conclusion and implications 

We have elucidated the origin of SLRs in the early solar sys- 
tem by developing a new model for 26 Al that relies on a physical 
mechanism (collect + injection, col lapse) similar in esse nce to 
the one presented recently for 60 Fe (iGounelle et al.ll2009l) . This 
new mechanism occurs naturally within a common mode of star 
formation, namely t hat of sequential star for mation within a gi- 
ant molecular cloud dHennebelle et al]l2009h . Within a few mil- 
ligrams of meteorites, SLRs therefore record physical mecha- 
nisms observed in the sky on scales varying from hundreds down 
to 1 pc. 

The identified sequence of events establishes the genealogy 
of the solar system. Our Sun is the great-grandson of a star com- 
plex (generation 0) containing 10s of thousands of stars, the 
grandson of a large GMC core (generation 1) containing a few 
thousand stars, and the son of a massive (>32 M©) star belong- 
ing to a cluster (1000-2000 stars with a preferred value of 1200 
stars, generation #2) born later within the same GM C. Assuming 
30 % star formation efficiency (lLada & Ladall2003l) and with an 
average stellar mass of M* =0.5 M©, our Sun was born to- 
gether with ~ six hundred fellow stars in its natal 1000 M© shell. 
Formation of the Sun in a relatively small clusteris in agreement 
with dynamical requirements, i.e. stability of planetary orbits, 
existence o f the Kuiper belt object Sedna, an d formation of the 
Oort cloud (lAdamsll2010l:lBrasser et al.ll2012l) . 

Given the size of its cluster at birth (six hundred stars), it 
is not inconceivable that our Sun was coeval to a massive star. 
Since the average number of massive stars in a cluster of size 
N is fsN x N, with fsN = 2.3 x 10" 3 , clusters with six hundred 
stars might contain between 1 or 2 massive stars, most probably 
B stars, leaving room fo r disk photo evaporation as suggested by 
iThroop & Bailvl (l2005h . 

It is expected that the few hundred stars born together with 
our Sun will share the same chemical and isotopic properties. 
These stars were true twins of our Sun Given that stellar clus- 
ters dissipate on timescales of 100 Myr (I Adams et al.ll200ll) . that 
revolution timescales around the Galactic center are on the order 
of 200 Myr, and that stars can r adially migrate over a few kpc in 
the Galaxy (Roskar et al. 2008), these fellow stars are now in to- 
tally unrelated places. Because the 26 Al enrichment mechanism 
we have identified is generic, many other stars are, however, ex- 
pected to contain 26 Al at a level close to that of our solar system. 
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